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Abstract 

We note that the new unitarity regime when scattering amplitude goes 
beyond the black disc limit (antishadowing) could help in the explanation of 
the regularities such as knee in the energy spectrum, existence of penetrating 
and long-flying particles and other features observed in the measurements of 
the extensive air showers which originate from cosmic particles interactions 
with the atmosphere. 
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The experimental and theoretical studies of cosmic rays are the important 
source of astrophysical information (cf. e.g. |.1|) and they simultaneously pro- 
vide a window to the future results of accelerator studies of hadron interaction 
mechanism at the LHC^ 
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Figure 1: Scaled energy spectrum of the cosmic rays, figure from Q. 

It can happen that the investigations of cosmic rays will give us a clue that 
the hadron interaction and mechanism of particle generation is changing in the 
region of ^/s = 3 — 6 TeVOHUl. Indeed, the energy spectrum which follows 
simple power-like law F(E) = cE~'^ changes its slope in this energy region and 
becomes steeper: index 7 increases from 2.7 to 3.1. It is important that the knee 
in the energy spectrum appears in the same energy region where the penetrating 
and long-flying particles also start to appear in the extensive air showers (EAS): 
the absorbtion length is also changing from A = 90 gjcvr? to A = 150 g/crn^ (cf. 
|l3l). There is also specific feature of the events at the energies beyond knee such 
as alignment cf. [i5j| and the references for the earlier papers therein. The above 
phenomena were interpreted as a result of appearance of the new particles which 
have a small inelastic cross-section and/or small inelasticity. These new particles 
can be associated with a manifestation of the supersymmetry, quark-gluon plasma 

'it should be noted that the value of the total cross-section extracted from cosmic rays measure- 
ments significantly depend on the particular model for elastic scattering, because measurements 
of the extensive air showers provide information on inelastic scattering cross-section onlvli2l. 
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formation and other new mechanisms. However, there is another possibility to 
treat those cosmic rays phenomena observed in EAS as the manifestations of the 
new unitarity regime (antishadow scattering mode) at such energies [6|. 

Unitarity of the scattering matrix SS^ = 1 implies, in principle, an existence 
at high energies s > sq, where sq is a threshold^ of the new scattering mode — an- 
tishadow one. It has been revealed in and described in some detail (cf. [8| and 
references therein) and the most important feature of this mode is the self-damping 
of the inelastic channels contributions at small values of impact parameter — an- 
tishadowing. The antishadowing leads to P(s, 6 = 0) ^ 1 at s — > oo, where P 
is a probability of the absence of the inelastic interactions, P(s, b) = \S{s, 6)p, 
where S is the elastic scattering S-matrix. 

Self-damping of the inelastic channels leads to asymptotically dominating role 
of elastic scattering. The cross-section of inelastic processes rises with energy as 
In s, while elastic and total cross-sections behave asymptotically as In^ s. The 
antishadow scattering mode could definitely be observed at the LHC energies and 
studies of the extensive air showers originated from the cosmic particles interac- 
tions with the atmosphere provide evidence for it as we will argue in what fol- 
lows. Starting at some threshold energy sq (where amplitude reaches the black 
disk limit at b = 0), antishadowing can occur at higher energies in the limited 
region of impact parameters b < R{s) (while at large impact parameters only 
shadow scattering mode can be realized). 

The inelastic overlap function r](s,b) becomes peripheral when energy goes 
beyond s = sq (Fig. 2). At such energies the inelastic overlap function reaches its 
maximum value atb = R{s), where R(s) is the interaction radius, while the elastic 
scattering occurs at smaller values of impact parameter, i.e. {b'^)ei < {b'^)inei- Note 
that 

1 f°° ^^rr. 



where i = tot, el, inel and 



T p/ i\ ^ da tot 2 ^ daei , , 1 daimi 

and unitarity condition in the impact parameter space is the following 

Im/M) = + r/(s,6), 

^Model estimates show that new scattering mode starts to develop right beyond Tevatron ener- 
gies, i.e. at ^/sq ~ 2 TeV [9J, which corresponds to the energy in the laboratory system £' ~ 2 
PeV. 
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where f(s, b) is the elastic scattering amplitude. The quantity (6^) is a measure 
of the reaction peripherality. Despite that the asymptotics for a^i and ainei are 
different, the quantities {b'^)ei and {b'^)inei have the same asymptotical energy de- 
pendence, proportional to In^ s. 

So, beyond the transition energy range there are two regions in impact pa- 
rameter space: the central region where self-damping of inelastic channels occurs 
(antishadow scattering at 6 < R{s)) and the peripheral region of shadow scattering 
at 6 > R{s). 




Figure 2: Impact parameter dependence of the inelastic overlap function in the 
framework of the unitarization scheme with antishadowing. Arrows indicate the 
directions of movement of minimum at 6 = and maximum at 6 = R{s) with the 
energy increase. In the region of 6 = R(s) the complete absorbtion takes place, 

i.e. \S{s,b = R{s))\^ = 0. 



At the energies s ^ sq small impact parameter scattering is almost elastic one. 

Thus head-on colliding particles will provide appearance of penetrating long- 
flying component in the EAS and such particles will spend only small part of 
their energy for the production of secondaries. The head-on collisions will lead 
to smaller number of secondary particle and it will provide faster decrease of 
the energy spectrum of cosmic rays, i.e. it will result in the appearance of the 
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knee. This qualitative picture will be explained in more detail in what follows. It 
should be noted that this effect has a threshold in the energy dependence. It is also 
important to note that due to small probability of the sequential head-on collisions 
the number of events with penetrating particles also should be small. Nontheless, 
such events have been observed in the experiments PAMIR ifTDIl . 

Antishadowing leads to suppression of particle production at small impact pa- 
rameters: 



i.e. multiplicity distribution 



and mean multiplicity n{s,b) in the impact parameter representation have no ab- 
sorptive corrections, but peripherality of darnel/ db'^ leads to suppression of parti- 
cle production at small impact parameters and the main contribution to the integral 
multiplicity n{s) comes from the region of 6 ~ R{s) (Eq. (HJl). This would lead 
to the events with alignment observed in EAS and also to the imbalance between 
orbital angular momentum in the initial and final states since particles in the final 
state will carry out large orbital angular momentum. To compensate this orbital 
momentum spins of secondary particles should become lined up, i.e. the spins of 
the produced particles should demonstrate significant correlations when the anti- 
shadow scattering mode appears ifTTl . Thus, the observed phenomena of align- 
ment in EAS 1^ and predicted spin correlations of final particles should have a 
common origin. The model estimate for the primary energy when these phenom- 
ena should appear is E'o — 2 PeV — Eq is the energy when the new unitarity 
regime starts to develop at small impact parameters. 

The detected particle composition of the EAS is closely related to the quantity 
known as gap survival probability. Antishadowing leads to the nonmonotonous 
energy dependence of this quantity (TT\. The gap survival probability, namely the 
probability to keep away inelastic interactions which can result in filling up by 
hadrons the large rapidity gaps, reaches its minimal values at the Tevatron highest 
energy and this is due to the fact that the scattering at this energy is very close to 
the black disk limit at 6 = (Fig. 3). It is clear that its higher value means higher 
fraction for diffractive component and consequently the increasing of this com- 
ponent would result in the enhancement of the relative fraction of protons in the 
observed cosmic rays spectrum. Otherwise, decreasing of this quantity will lead to 




(1) 
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Figure 3: Energy dependence of gap survival probability. 

increase of pionization component and consequently to the increasing number of 
muons observed as multi-muon events. Experiment reveals that relative fraction 
of protons in cosmic rays also shows nonmonotonous energy dependence (cf. Fig. 
4). To explain such dependence an additional component is introduced ad hoc at 
the energies above 3 ■ 10'' GeV. It was shown that account of the antishadowing 
makes an introduction of this ad hoc component unnecessary. 

The inelasticity parameter K, which is defined as ratio of the energy going 
to inelastic processes to the total energy, is important for the interpretation of the 
FAS cascades developments. Its energy dependence is not clear and number of 
models predict the decreasing energy dependence while other models insist on 
the increasing energy behaviour at high energies (cf. e.g. ifT^ ). Adopting simple 
ansatz of geometrical models where parameter of inelasticity is related to inelastic 
overlap function we can use the following equation for {K) lfT4ll 

to get a qualitative knowledge on the inelasticity energy dependence. The esti- 
mation of inelasticity based on the particular model with antishadowing [9| leads 
to increasing dependence of inelasticity with energy till ~ 4 ■ 10^ GeV. In this 
region inelasticity reaches maximum value (K) = 1, since aei/crtot = 1/2 and 
then starts to decrease at the energies where this ratio goes beyond the black disk 
limit 1/2. Such qualitative nonmonotonous energy dependence of inelasticity is 
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the result of transition to the antishadowing scattering regime. The distribution 
on the inelasticity is related to the distribution on the effective mass number, i.e. 
changes of A are equivalent to changes of (K), and, for example, high-inelasticity 
primary proton interaction produces the same result at the ground level as the low- 
inelasticity primary interaction of the heavy nuclei ifT?! . The available experimen- 
tal data on the average logarithm of the effective nuclear mass number, extracted 
from the energy dependence of the depth of EAS maximum, have large error bars, 
but they also indicate a nonmonotonous energy dependence with the maximum in 
the region ^ ~ (4 - 5) ■ 10^ GeV fl\. 

It is also worth to note that the maximum in inelasticity energy dependence, 
when the pionization component is maximal, is correlated with the minimum of 
the relative component of protons in the EAS, the following simple relation can 
be supposed 

i.e. the relative proton component in the detected EAS should have a non-monotonic 
energy dependence and this is in agreement with the experimental analysis repre- 
sented in Fig. 4. 
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Figure 4: Relative fraction of protons in EAS, figure is taken from 0. 

It should be noted that the behaviour of the ratio cfei/crtot when it goes to 
unity at s ^ oo does not imply decreasing energy dependence of ainei- The 
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inelastic cross-section ainei increases monotonically and it grows as In s at s — > 
oo. Such a dependence of aimi is in good agreement with the experimental data 
and, in particular, with the observed falling slope of the depth of shower maximum 
distribution [.16J . The predicted numerical value of the inelastic cross-section is 
ainei{s) ~ 76 mh at the LHC energy -^s = 14 TeV. This value is also in a good 
agreement with the value for this quantity extracted from the proton-air inelastic 
cross-section |7|. This approach provides a reasonable description lITTll of the 
energy dependence of mean multiplicity and leads to its power-like growth with a 
small exponent. 

The relation of the knee and other effects observed in the EAS measurements 
with the modification of particle generation mechanism is under discussion since 
the time when they were discovered. We propose here one particular realization of 
this idea — an approach where the corresponding particle generation mechanism 
in EAS is strongly affected by the unitarity effects and the energy region between 
the knee and the ankle coincides with the transition region to the scattering mode 
where antishadowing develops at small and then at moderate values of impact 
parameter, i.e. the energy spectrum of the primary cosmic particles Fq{E) is 
modulated by the significant variation of the scattering matrix 5* in the energy 
region starting from about Ei ~ 10^ GeV and finishing at about E2 — 10^ GeV 
and this resulting in the regularities in the observed spectrum F{E) measured 
in the EAS studies. Below the energy Ei and beyond the energy E2 variation 
of scattering matrix is slow and the primary energy spectrum Fq is almost not 
affected. It seems to be a rather natural explanation of the observed regularities 
in the EAS measurements and has a close interrelation with the nonmonotonous 
energy dependence of gap survival probability and inelasticity. This hypothesis is 
based on the saturation of the unitarity and can be experimentally checked at the 
LHC [8|. The studies of the proton scattering in the forward region at the LHC 
will be very helpful for improving the interpretation of the results of the cosmic 
rays experiments. 
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